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The thermodynamical properties of a pure Yang-Mills theory SU (N) is described by a gas of glu- 
onic quasiparticles with temperature dependent mass m(T) and a bag function B(T). The analytic 
behavior of m(T) and the pressure p in the temperature range 2.5-5X' c are derived and constraints 
on the parameters defining B(T) are discussed. The trace anomaly 8 = p — 3p is evaluated in the 
high T domain: it is dominated by a quadratic behavior 9 — nKT 2 , where n = 2(N 2 — 1) is the 
number of degrees of freedom and K is an integration constant which does not depend on the bag 
function B(T). The quadratic rise of 6 is in good agreement with recent lattice simulations. 



I. INTRODUCTION AND SUMMARY 



The concept of quasiparticle is a valuable tool toward an effective description of complicated interactions. An 
important example is that of Yang-Mills (YM) SU(N) theories at nonzero temperature T, where their intrinsic 
nonperturbative nature renders the applicability of a perturbative treatment limited [l|. Descriptions of the YM 
system at T > via gluonic quasiparticles have been performed in a variety of approaches, e.g. Refs. 041] and 
refs. therein. Different Ansatze have been tested and the outcoming thermodynamical quantities, such as energy 
density and pressure, have been compared to lattice results [9l4Tll]. An appropriate, but at the same time analytically 
tractable, description of a gas of gluons is also a necessary step toward the understanding of the quark gluon plasma, 
see for instance Ref. [l3[ for a review. 

Besides the temperature dependent gluonic mass m = m(T), also a temperature dependent bag energy B = B(T) 
has been widely used to describe nonperturbative properties of YM theories, such as the trace anomaly and the gluon 
condensate. With these two basic ingredients the energy density and the pressure read (see Refs. [3, 3 and refs. 
therein) : 



p = Pp + B(T) ,p = p p -B(T) 
where the suffix 'p' denotes the quasiparticle part: 
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where J k = J = Jo°° ITtt^ ana - n re P resen ts the degeneracy of the particle spectrum. Restricting to perturbative 
degrees of freedom only, one has n = 2 (TV 2 — 1) in the case of a SU{N) YM theory. 

In this article we study the high T properties of the gluonic gas expressed in Eqs. ([1} and ([2]). On the practical 
side, we are interested in the behavior in the range between -say- 2.5T C and 5T C , where T c is the critical temperature 
for deconfinement, above which gluonic quasiparticles are the relevant degrees of freedom. In this energy range lattice 
data for the energy density and pressure show a plateau (9l-[il|. which is compatible with a linear increase of the 
quasiparticle mass with the temperature T, m oc T. Because of this linear increase the Stefan-Boltzmann limit is not 
reached, but a saturation at lower values is obtained. This situation persists as long as logarithmic corrections due 
to the running coupling of QCD are neglected. In fact, the correct perturbative behavior is m oc T j y/log(T j Aqcd), 
which implies a slow approach of the pressure and energy density to the expected Stefan-Boltzmann values, see the 
lattice simulation in Ref. [12J. In the present work we do not include these logarithmic corrections. This simplification 
is applicable in the temperature range 2.5-5T c , which is high when compared to the critical temperature, but not high 
enough for the logarithmic terms to become relevant. 

Strong constraints on the gas of quasiparticles can be derived by imposing that the system fulfills thermodynamical 
self-consistency 04^] , which is a consequence of the first principle of thermodynamics: 



P = T 



dp 
dT 



The bag constant B(T) is assumed to have the following behavior: 

B(T)=B NP (T) + B P {T) 



(3) 



(4) 
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with B NP {T) = ncT a and B P (T) = nST 4 for 2.5T C < T < 5T C . (5) 

The term Bp(T) = nST 4 is included in order to recover the expected results of perturbation theory in the high 
temperature regime, where the effective gluon mass grows linearly with T (up to the previously mentioned logarithmic 
corrections, which are not considered here.) 

The term B^piT) = ncT a , where a is a real number smaller than 4 and c is a constant with the dimension of 
[Energy] 4- ", describes the 'nonperturbative' bag contribution relevant above the phase transition [IJj]. At the present 
stage the only and general requirement about the function B NP (T) is that it is dominated by a power-like term T a in 
the high T region. It is the aim of this work to constrain the value of a and c by using mathematical considerations 
and lattice results. 

As a last remark we stress that the bag function B(T) = Bnp(T) + Bp(T) is proportional to the number of degrees 
of freedom n — 2(N 2 — 1), in agreement with general large N scaling arguments [l 51 ] . 

The temperature dependent mass m(T) can be analytically evaluated at high T (details are in Sec. 2) and takes 
the form (a ^ 2): 



/ 4tt 2 ac A 2 8tt 2 (5 

m (T) = TJ 7 a ° T<*-4 + k^- ^-L , (6) 
y D(a ) 2 - a T z D(a ) 

where -D( a o) is a positive real number which shall be specified later on, a, c, S are the already introduced parameters 
defining B(T), n is the number of degrees of freedom, A is the Yang-Mills scale and finally k is an integration constant 
related to the differential equation ([3]l: k is not determined by the choice of B(T) but is a further parameter entering 
in the model. 

Mathematically and physically based considerations about Eq. ([6]) will lead us to establish nontrivial relations 
between the parameters: 

• The parameter 8 parametrizes the linear increase of m with T. As evident from Eq. (|6]), 5 < in order to avoid 
an imaginary mass at high T . Moreover, a careful study of the equations will lead us to establish also an upper limit 
on its absolute value: \8\ < 0.0151. 

• Many approaches show that the nonperturbative bag function B^p is a positive number [l5|, [Hj]. In the present 
framework it means that c > 0. Moreover, the contribution of nonperturbative physics to the effective gluon mass 
should be positive, a property which also assures that no instability emerges at low T. Thus, it is a general physical 
requirement that the nonperturbative contribution to m 2 (T) given by T^)^T a - 2 in Eq. © is also positive. 
This, in turn, implies a consistent limitation on the choice of the parameter a: < a < 2. Note, for the very same 
requirement we are led to conclude that the integration constant k is positive. 

• The case a ~ 2 is somewhat peculiar because of the emergence of logarithms in the solution, see details in Sec. 
2.4. However, it is also unfavoured because of similar arguments. 

As a next step of this work we turn to the explicit expression for the pressure p in the high T limit, see Sec. 3 for 
details. We shall find that p is expressed by the sum of three terms (a ^ 2): 

2 n 

p = - nc ^—^ Ta - ^ D M kA2T2 + (fpM - nS ) Ti ■ (7) 

The first, negative term scale as T a (just as Bjyp); the second, also negative, term scales as T 2 (in agreement with the 
phcnomenological argument of Ref. [13]) and is proportional to the constant k\ the third, positive term goes as T 4 , 
but the coefficient is slightly smaller than the Stefan-Boltzmann limit, see later on for details. The overall pressure 
scale as n oc N 2 in agreement with large N scaling arguments. 

The final and main subject of the present work is the study of the trace anomaly at high T (Sec. 4). The trace 
anomaly 9, defined as 

= p - 3p = AB + p p - 3 Pp , (8) 



has been in the center of a vivid debate in the last years (see Refs. |17H22| and refs. therein). Strict dilatation 
invariance would imply that 9 vanishes in a dilatation invariant theory, such as a gas of photons. In a YM theory this 
symmetry is broken by quantum effects and 9 does not vanish: this is the so-called trace anomaly. We aim to show 
that, in the context of a gas of quasiparticle with the general form of the bag constant B(T) given in Eq. (0}, the 
following high T behavior holds: 

9 = p-3p = nCT a + nKT 2 for T > 2T C , (9) 

where C = 2c|5§ is a constant determined by the 'nonperturbative' parameters of the model a,c (i.e., those pa- 
rameters which define the nonperturbative bag function Bnp(T) — ncT a ), and K — fcA D ( a °) j s proportional to the 
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previously introduced integration constant k. The main result is that the trace anomaly 9 can be decomposed in a 
term which behaves as the 'nonperturbative' contribution to the bag constant B^p{T) — ncT a , and a term which 
goes as T 2 . Restricting to the favoured interval < a < 2, one is led to conclude that the quadratic rise dominates 
at T large enough. Remarkably, the T 2 rise of 9 is a general property, which is independent on the nonperturbative 
bag constant -Bjvp. 

It is indeed remarkable that a quadratic rise of the trace anomaly, 9 ~ aT 2 , is found in Ref. [17| . where an analysis 
of the lattice data of Ref. [Toj has been performed. Later on, this quadratic rise has been confirmed in recent lattice 
works [ll|, [23|. In particular, in the lattice study of Ref. |23| the trace anomaly has been investigated for various 
pure Yang-Mills theories SU(N), N = 3, 8. The behavior 9 oc nT 2 (with a direct proportionality to the degeneracy 
number n = 2(N 2 — 1)) is indeed found for N = 3, 8 in a range between 2 and 5T C . Thus, the result of Eq. © 
may explain in a natural way the emergence of such quadratic behavior of 9 at high T . 

In the following sections we present the detailed derivations of the outlined results: in Sec. 2 and 3 we derive the 
expressions for m(T) in Eq. (|6]) and for the pressure p in Eq. ([7]). In Sec. 4 we present the calculation leading to the 
expression of the trace anomaly 9{T) in Eq. ([9]). Finally, in Sec. 5, we briefly outline our conclusions. 



II. TEMPERATURE-DEPENDENT QUASIPARTICLE MASS m(T) 



A. Differential equation for m(T) 



In order to obtain the differential equation for the quasiparticle mass m = m{T), we plug the expressions of Eqs. 

dp 
dT 

1 1 



(HJ) and © into the the thermodynamical self-consistency relation p — T^ — p 



dB dm 2 

— —nl(m) — , l(m) = 
dT y ' dT ' y ' 



where the bag function B(T) is given in Eq. Qj. 
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B. Use of dimensionless functions 



It is convenient to rewrite the equations by using dimensionless quantities. To this end we introduce the dimen- 
sionless temperature 

an 

where A is the Yang-Mills scale, which is of the same order of the critical temperature T Cl A ~ T c . 
The dimensionless particle contribution to the energy density and pressure 

Pp = -fl . fp = 7fi ( 12 ) 

read 



P » = ^L ^^ri' (13) 



^ = -^l dx(x 2 ln(l-e-^)) . (14) 
The function a = a(\) is the 'dimensionless mass': 

a ^ a{x) ^ain^!^m. (15) 

T AA v ' 

We also define the dimensionless constant 7 as 

7 = ck a ~ A . (16) 
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In this way the ground-state dimensionless energy density p gs and pressure p gs 



= _ = B(T) 

Pgs Pgs rp± 

read explicitly (see Eq. ((4]) 



Pgs — Pgs — rp± — rp4 \ L ' ) 



P gs = -p gs =n 1 \ a - 4 + nS . (18) 
The full dimensionless energy-density and pressure read 

P=Pp+Pgs -.P^Pp+Pgs ■ (19) 

The thermodynamical self-consistency of Eq. ([3]) can be rewritten in terms of the reduced energy density and pressure 
as 

P-3 P =a|. (20) 
In terms of the dimensionless quantities he differential Eq. (|10[) takes the form: 



xa -4^ A A n d(a?\ 2 ) D(a) 
na 7 A +4n<5=- — — — , (21) 



2tt 2 dA 2A ' 

where the integral D(a) is given by 



D(a) = / dx ■ -= = . (22) 

J a Vx 2 + a 2 e V^+^ _ i v ; 

Note, the dependence on the degeneracy number n = 2(N 2 — 1) factorizes, so that the equation for a(A) is independent 
on the number of colors N. This is in agreement with the general expectation of large N scaling, according to which 
the effective gluon mass m(T) scales as N°. 



C. Constraints on the parameter S 

In the limit of large A (i.e., large T) one has a(A >> A c ) — > do. In this way, besides the logarithmic corrections, the 
effective mass exhibits a linear growth m = anT, in agreement with the expectation of perturbation theory [24| and 
with high T effective approaches, e.g. Ref. [25| . In the present phenomenological approach the numerical value of ao 
is related to the parameter 5 by studying the asymptotic behavior of Eq. (|2T|) : 

8 = ~<%D(ao). (23) 

In Fig. 1, left panel, the quantity S is plotted as function of oq. The two properties mentioned in the Introduction 
can be easily proven: S < and \5\ < max (^an-P(an)) = 0.0151. 



Various lattice simulations of Yang-Mills system [lCj, 111] show that the Stefan-Boltzmann limit of the energy density 
and the pressure 

Psb = n 7^ , Psb = n— (24) 

is not reached at 5T C . On the contrary, a saturation at a lower value of about 90% of the Stefan-Boltzmann limit is 
observed. Such a saturation is obtained in the present quasiparticle approach by a nonzero value of ao (i.e., a nonzero 
value of S). For high A the function p(A)/p SB approaches the asymptotic value (p p (a ) — nS)/p SB (see Sec. III). In 
Fig. 1, right panel, the quantity (p(X)/p SB ) x>>1 — (p(A)/psb)a>>i is plotted as a function of ao. In order that at 
high A the ratio p(A)/p SB ~ 0.9 holds, the value an — 0.83 is required. This, in turn, implies that 5 ~ —0.0070. Note, 
similar values for ao have been obtained in Refs. [HQ. 

As explained in the Introduction, the Stefan-Boltzmann limit for the energy density and the pressure is reached at 
much higher temperatures [l2j |, at which the logarithmic decrease of a(A) becomes relevant. 
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FIG. 1: Left panel: the parameter S is plotted as function of the asymptotic value a(A — » oo) = ao- Right panel: the ratio 
(p/Pss)a>>i = (p/Pss)a>>i is plotted as function of ao. The dot corresponds to p/p SB = 0.9 and ao = 0.83; this is the 
saturation value obtained in lattice simulations 1231 ■ 
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FIG. 2: Comparison of the analytic, approximate expression of Eq. (|25[) (solid line) with the full solution of the differential Eq. 
(f2"T) (dashed line) for a particular numerical case: a = 1, 7 = 0.034, k — 10. 



D. Analytical solution a(A) in the large A domain (i.e., m(T) in the large T domain) 



An analytical solution of Eq. (|2ip can be obtained in the limit of large A by approximating the function D(a) by 
its asymptotic values D(ao). In this limit Eq. (|2"Tj) can be easily solved and one obtains for a 2 (X) (a 7^ 2): 



In Fig. 2 the analytical expression (|2"5|) is compared to the numerical result of the differential Eq. (|2"Tj) for a 
particular numerical choice; the numerical solution takes into account the full dependence on the function D{a) 1 
which is ignored in the analytical one. However, the qualitative agreement, which becomes better and better by 
increasing A, is visible. 

The analytical expression of the function m 2 (T) is easily obtained by making use of Eqs. (ITT|) . (fT5)) and (|25l) : 
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As a result of the obtained expression for m(T) we can discuss the constraints on the parameters c and a: 

(i) The squared mass constitutes of three terms. The first term in Eq. (|26[) depends on the parameters c and a, 
which define the nonperturbative bag pressure B^p — cT a . The second contribution in Eq. (|26|) . fcA 2 , is constant 
and is directly proportional to the integration constant k. The last term in Eq. (|26p dominates at high T, implying 
the linear behavior m(T) ~ a$T ', as predicted by pcrturbative calculations and effective theories of QCD. As already 
anticipated in Sec. 1, the natural requirements that Bnp — cT a and the corresponding contribution to m 2 (T), given 

D(ao) T^aT a ~ 2 , are positive numbers implies that: 7 > 0, < a < 2. Similarly, k > 0. 

(ii) In the limit a = the first term in Eq. (|26|) vanishes and the second, constant contribution dominates: this 
situation corresponds to the simple approximation with a constant gluon mass. This possibility has however been 
ruled out by a precise comparison with lattice data More generally, when a is small, the first term in Eq. (|26[) 
is also negligible (unless the parameter c is anomalously large); at low T only the second constant term survives. On 
the contrary, when a ~ 2 the first term is very large, unless the parameter c is very small. We thus conclude that a 
should be not to close to the boundaries and 2, but somewhere in between. 

(iii) In the case a — 2 & slightly different solution is obtained: 



in 2 -f log A k 2 
D(ao) A 2 +A 2 



« 2 (A) = -7^7TT-^r + T2+ a o, ( 2 7) 



where an additional logarithm arises. The mass contribution of the nonperturbative first term is negative for 7 > 
(that is for the here considered choice c > 0). Thus, also the case a = 2 is regarded as phenomenologically unfavoured. 

(iv) In order to include the logarithmic corrections in the very high T domain one should modify the perturbative 
bag pressure Bp(T,) in such a way that the perturbative mass behavior m cx T/ylog T/A holds. One obtains the 
constrain T~ 3 dBp/dT cx (21n _1 A — 21n~ 2 A) , thus leading to more complicated expressions involving logarithms. A 
detailed study of this subject represent an interesting outlook. Although the formulas will be more involved, a link 
with studies of Ref. Q, in which the starting point is the perturbative behavior at very large temperature, can be 
driven. 



III. PRESSURE IN THE HIGH T DOMAIN 



We turn to the explicit expression of the pressure for large T. To this end we expand Eq. (|I4[) around the asymptotic 



"o- 



/ dp A a) \ 

p p (a)=p p (a )+{-^) (o-oo) + ... (28) 
Using the equality ( J = — 2^2 a o-C( a o) an d approximating ao(a — a ) = a 2 — a 2 , (valid at the considered order) 

V a J ao 

one gets 



The full dimensionless pressure p = p v + p gs reads at high A: 



p = -n 7 -^-A- 4 - ^D(ao)^ + {p p (a ) - nS) . (30) 



2-Q 47T 2 A 

By multiplying by T 4 we find the pressure p for large T: 

p= - nc 2h^ Ta ~£^ D{ao)kA2T2 +(p P ( a o)- nS ) T4 - ( 31 ) 

We thus have also decomposed the pressure into three contributions: the first term in Eq. (l3"Tj) scales exactly as the 
bag function B^p- The second, negative term in Eq. (1311) scales as T 2 and is proportional to the integration constant 
k. Note, a similar quadratic contribution to the pres"re has been postulated In Ref. 0. There is, however, an 
important point to stress: here we have shown that there is is no need to introduce at hand a quadratic contribution 
to the pressure from the very beginning. The quadratic contribution naturally emerges as the result of the equation, 
independently on the choice of the bag function Bnp(T). Finally, the last term in Eq. (|31[) describes the high T 
asymptotic limit, which differs from the Stefan-Boltzmann value as depicted in the right panel of Fig. 1. For other 
works on the pressure in the high T domain see Refs. [26l . [27| and refs. therein. For a direct comparison with lattice 
data see Fig. 3. 
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IV. TRACE ANOMALY 



The interaction measure 



A = — = p-3p (32) 
is evaluated by making use of the thermodynamical self-consistency of Eq. (|2"U|) : 

A = p - 3p = A— = A— + j (33) 

where in the last step the dimensionless pressure has been decomposed into its particle and ground-state contributions. 
The ground-state contribution is easily evaluated: 

A -jf =«7(4-«)A Q - 4 . (34) 
d\ 

The calculation of the particle contribution to the interaction measure goes via two steps. First, we rewrite it by 
making use of Eqs. (|T4"|): 

, dpr, n . „ , . da n , „ , . da 2 

x ^ = -^ XD ^ a dx = -^ XD ^- ™ 

As a second step, in the large A domain one can replace D(a) with the asymptotic value D(ao) and then evaluates 
the derivative by using Eq. (j2"5)l : 

dp 4- a 4 L>(a ) fe 

A dA = " a7 2^ A +n ^A^- (36) 

By putting the results together one finds for A > 2A C (a ^ 2) 

A = p-3p = 2n 7 ^A- + n^' , (37) 
2 — a 27r z A z 

which consists of two terms: a term which scales as the ground-state contribution p gs , and a term which scales as 
A- 2 . 

By multiplying A by T one obtains the trace anomaly 9 as function of T (a ^ 2): 

6 = 2nc—T a + n^^kk 2 T 2 , (38) 
2 — a 2n 2 

which is Eq. ^ discussed in the Introduction: the first term depends on c and a (i.e., the parameters which define 
Bnp{T)), while the second term describes a quadratic rise of 9, is proportional to the integration constant k and is 
independent on the bag function B. The very same term proportional to k was responsible for a constant contribution 
to the effective gluon mass, see Eq. (|2l>|) . In the favoured range < a < 2 the rise T a is realized for small T, while the 
quadratic rise of 9 dominates for large enough temperature. The temperature at which this change happens depends 
on the particular numerical values of the parameters, and cannot be determined by analytical considerations. 

It is however possible to use some lattice results about the trace anomaly in order to constrain the numerical values 
of the parameters of the model. The quantity 9 scales as 3.3T C 2 /T 2 for T > 1.5-2T C 4]. Then, from Eq. and 
A ~ T c it follows that k ~ 10. It is also possible to obtain a rough estimation of the upper limit of the parameter 
7 = cA 4 ~" ~ cT^~ a - In fact, the quantity 9/T 2 is, to a good approximation, constant for T > 2T C [17|. This, in 
turn, means that the first term in Eq. (|3"8")l is smaller than the second term for T > 2T C . Through simple algebra one 
obtains the upper limit 7 < 2 ~_ a a D ^ k 2 1 - a . 

We now turn to a direct comparison of our theoretical curves with the lattice results of Ref. Since our 

theoretical functions depend on the variable A = T/A, where A ~ T c but not exactly equal, care is needed: it is 
first necessary to determine A. To this end we chose A in such a way that the theoretical result for the interaction 
measure A reproduces the lattice point at the highest simulated value of T/T c = 4.57 (at which A = 0.10). One 
obtains the relation A = 1.55T" C . In Fig. 3 the plot of the interaction measure (left panel) and energy and pressure 
(right panel) are shown: it is visible that the agreement is acceptable for T > 2.5T C and increases for increasing T. 
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(a) Interaction measure (b) Pressure and energy density 

FIG. 3: Comparison with the lattice results of Ref. The parameters a = 1, 7 = 0.034, k — 10 are used. Left panel: the 

interaction measure is plotted. The lattice point for T/T c = 4.58 has been uused to determine A = 1.66T C . Right panel: the 
quantities p (the upper curve is the theoretical result, the dots the lattice points) and 3p (the lower curve is the theoretical 
result, the triangles are the lattice points) are shown. 



On the contrary the theoretical results for T/T c < 2.5 deviate from the lattice simulations. This is expected because 
the present version of the model cannot describe the physical properties close to the phase transition. 

Further comments are in order: 

(i) For the here considered case c > (i.e. for B^p > 0) a positive contribution of the first term to the trace 
anomaly is obtained for < a < 2 (see Eq. (|38|l). which is the same interval outlined previously. If, on the contrary, 
2 < a < 4 one would have a negative 9 for high enough T, in disagreement with all present simulations. This 
represents a further confirmation of the outlined range of a. 

(ii) The case a = 2 leads to a slightly modified form: 



e = 4nT'[c-clog^jj+-^kA 2 j . (39) 

A negative 9 at high T is realized. This fact is at first sight a further argument against the choice a = 2. However, 
the appearance of the logarithms in the solution implies that a full study of the present case is only possible when 
the logarithmic corrections are taken into account. 

(iii) The behavior of 9 in the large T domain as measured on the lattice is still subject to an ongoing discussion. 
In the work of Ref. [IH , also based on the lattice data of Ref. [Io| , it is found that 9 growths linearly (rather than 
quadratically) with T : 9 = aT for 2.5T C < T < 5Tc, where a ~ 1.5 GeV 3 in the SU(2)case and a ~ 1.7 GeV 3 in the 
SU(3) case. We also refer to the similar results obtained in the older works of Refs. [2^, [29|. In Ref. [3(| the linear 
rise has been confirmed by studying the lattice data of Ref. [ll| . Our result ([9]) can indeed also account for an initial 
non-quadratic behavior of 9 (the linear one being realized for a = 1, see [l9|), which persists as long as the quadratic 
term does not become dominant. 



(iv) The linear rise of 9 has been derived within the theoretical framework described in Refs. [31], and further 
investigated in Refs. [lS, 32]. The non-perturbative sector of SU (2) or SU(3) YM theories is described by a composite, 
(adjoint-)scalar field <p in the deconfined phase (T > T c ), which emerges as an 'average' over calorons and anticalorons 
(topological objects which correspond to instantons at nonzero T [33[) with trivial holonomy, see [3l| for a microscopic 
derivation and 32] for a macroscopic one. On a length scale I > \4>\ 1 it is thermodynamically exhaustive to consider 
only the average field 4> and neglect the (unsolvable) microscopic dynamics of all YM-field configurations, such as 
calorons and monopoles. One can then build up an effective theory for YM-thermodynamics valid for T > T c , in 
which the scalar field <f> acts as background field coupled to the residual, perturbative gluons. On a phenomenological 
level it contributes to the energy and pressure as a temperature-dependent bag constant jBatp = 47rA 3 T, i.e. with 
the parameters a — 1 and 7 = 4-7T. Note, in the theoretical framework of Ref. [31| the constant k was set to be very 
small, thus the quadratic rise starts to dominate only at very high temperatures and for this reason does not affect 
the phenomenology between 2.5T C and 5T C . The linear growth with T of the stress-energy tensor in the pure SU(2) 
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YM theory is obtained as [Til: 9 = p - 3p T - Tc 24ttA 3 T ~ (1.7 GeV 3 )T. The coefficient 1.7 GeV 3 is similar to 1.5 
GeV 3 found in Ref. [l8[. A similar result holds in the SU(3) case. 

(v) A linear growth of the trace anomaly 9 has also been obtained within the theoretical approach described in Ref. 

20], in which a Gribov-type dispersion relation is used. On the contrary, a quadratic rise of 6 is the result of Ref. 

21| . in which a dimension- two gluonic condensate is studied. 

(vi) In the present work we concentrated on the high T side. If we assume that a power-like behavior Bm p (T) — cT a 
is valid for lower temperature, the quasi-particle gluon mass blows up at a critical temperature T c (which depends 
on the numerical values of c and a, e.g. in Ref. |3l|). This fact may signalize a confinement /deconfinemcnt phase 
transition. For T < T c a quasi-particle description is no longer possible and the system should be described by 
different degrees of freedom (such as glueballs [34[ and other nonperturbative states). 



V. CONCLUSIONS 



In this work we have performed an analytical study of the high temperature properties of a gas of gluonic quasi- 
particles with a temperature-dependent bag function. The expression of the quasiparticle mass m(T), pressure and 
trace anomaly have been derived analytically for large T. 

The implications and constraints on the parameters of the bag function B(T) = B^p{T) + Bp(T) have been 
discussed: for the nonperturbative contribution Bnp(T) = ncT a (with c > 0) we have found the following constraint 
on the parameter a: < a < 2. This result follows from the requirements that the gluon mass m(T) does not become 
imaginary for decreasing T and that the sign of the trace anomaly 9 = p — 3p at large T is positive, in agreement with 
lattice simulations. The behavior of the trace anomaly 9 = p — 3p at high temperatures consists of two contributions, 
9 = n#T a + n#T 2 , with a first term which goes as T a , just as the the nonperturbative bag function, and a second 
term which goes as T 2 and dominates the high T behavior. It is remarkable that this quadratic contribution is general 
and does not depend on the choice of the bag function. Such a quadratic behavior, and also the expected scaling with 
the degeneracy number n, have been confirmed in the lattice simulations of Refs. [TT1, [23j| . 

A variety of improvements of the present approach represents an outlook for the future: (i) Inclusion of the 
logarithmic corrections in order to make contact with the very high T behavior, in which the pressure and the energy 
density slowly approach their Stefan-Boltzmann limit, (ii) On the low-temperature side a fit should be performed 
in order to determine the behavior of the nonperturbative bag pressure Bnp(T) beyond the simple power- law used 
in this work for analytical considerations, (iii) Calculation of viscosities following Ref. [35| can be performed, (iv) 
Inclusion of further degrees of freedom: quarks for T > T c and confined states (glueball and mesons) for T < T c . 
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